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Abstract 

The muon anomalous magnetic moment is investigated in the MSSM 
for tan/3 —)• oo. This is an attractive example of radiative muon mass gen¬ 
eration with completely different qualitative parameter dependence com¬ 
pared to the MSSM with the usual, finite tan/3. The observed, positive 
difference between the experimental and Standard Model values can only 
be explained if there are mass splittings, such that bino contributions domi¬ 
nate over wino ones. The two most promising cases are characterized either 
by large Higgsino mass jj, or by large left-handed smuon mass rriL. The 
required mass splittings and the resulting are studied in detail. It 

is shown that the current discrepancy in can be explained even in cases 
where all SUSY masses are at the TeV scale. The paper also presents use¬ 
ful analytical formulas, approximations for limiting cases, and benchmark 
points. 
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1 Introduction 


The muon anomalous magnetic moment = {g — 2)^/2 provides a tantalizing 
hint for physics beyond the Standard Model (SM). The current discrepancy 
between experiment and theor}|^ is 


exp 


,SM 


= (28.7 ±8.0) X 10 


-10 


( 1 ) 


It can be explained by a variety of new physics models below or at the TeV 
scale. Generally, the new physics contributions are strongly correlated with 
the loop contributions to the muon mass, , and they are suppressed by two 
powers of the typical new physics mass scale M^p (5, 19 
written as 


This relation can be 


NP _ ^ 

^^NP 


( 2 ) 


where Gnp = /rUfj,) is given by the model-dependent relative contribution 

to the muon mass. 

Of special interest are, therefore, models of radiative muon mass generation: in 
these models, amounts to the entire physical muon mass and Gnp = 0{1). 


NP 


Excluding hne-tuning in the muon mass, these models yield the largest , 
compared to other ones with the same new physics scale Mnp. As the estimate 
^NP _ holds, the observed deviation can in principle be explained 

by a new physics scale of the order 2 TeV. 

The general idea of radiative muon mass generation can be realized in the 
renormalizable and calculable framework of the minimal supersymmetric stan¬ 
dard model (MSSM). As discussed below in Sec. there are two distinct possi¬ 
bilities, one of which has already been studied in Refs, 
considered in the present paper is the limit 


20 22 . The possibility 


tan /? = — _). oo 
Vd 


( 3 ) 


with the up- and down-type Higgs vacuum expectation values v^^d- In this limit 
Vd, the muon mass, and all other down-type lepton and quark masses vanish at 
tree level. The masses arise from finite loop diagrams generating non-holomorphic 
couplings of down-type fermions to the “wrong” Higgs doublet Hu- These loop 
diagrams are also important for hnite tan/?, and they have been discussed ex¬ 
tensively in the literature, often in the context of R-physics 23 , but also in the 


^ The experimental result has been obtained at BNL |^; a fourfold improvement in pre¬ 
cision is expected from the new experiments [^|^. For reviews of the theory prediction see 
Refs. [^1^; recent theory progress has been achieved on the QED [6j|^, electroweak |^, and 
hadronic contributions 10 -14 ; for specific reviews of the hadronic light-by-light contributions 
and expected future improvements see 15 -18 . The value quoted in the text is based on the 
hadronic contributions from Ref. 
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context of the muon magnetic moment 24 . Because of these loop diagrams, the 


limit tan /3 —>■ oo exists and is phenomenologically viable 25,26 


We consider in the MSSM for tan/3 —>■ oo for two reasons. On the one 
hand, this scenario is a calculable realization of the generally interesting idea of 
radiative muon mass generation, and the results can be indicative of the more 
general situation. In particular, it is no special case of the model-independent 
analysis of Ref. 27 , where simplihed models with only two relevant particle 
masses have been considered. We will see that all cases of interest here involve 
at least three relevant masses at the TeV scale. 

On the other hand, the tan/3 —)• oo limit opens up an intriguing area of 
supersymmetry (SUSY) parameter space where behaves qualitatively very 
differently from the standard MSSM case (with moderate tan/3). This standard 


MSSM case is well known (for reviews see [28fp30] ; recent works are 31-4^), 
and it requires SUSY masses in the few-hundred GeV range in order to explain 
the deviation ([^; the LHC experiments, however, start to exclude parts of the 
relevant parameter space 31,46 . Thus it is well motivated to ask whether SUSY 


can explain the deviation even if all SUSY masses are much higher, at the TeV 
scale. The answer can be expected to be provided by SUSY radiative muon mass 
scenarios such as the one considered here. 

Our approach is a low-energy phenomenological one, with the aim to answer 
whether and how TeV-scale SUSY masses are compatible with the deviation ([^ 
in the MSSM limit tan/3 —)■ oo, also taking into account other experimental data 
and the internal consistency of the theory. Nevertheless it is important to briefly 
take a top-down perspective and review how inhnite tan/3 might result from a 
more fundamental theory. One appealing possibility is an unbroken continuous 
R-symmetry, which forbids the Bfi soft breaking term. Such a symmetry has 
been used in Ref. 47 to construct a “one Higgs doublet model”, in which Vd = 0. 


R-symmetric models such as this one however contain non-MSSM degrees of 
freedom, and therefore our MSSM study does not directly apply to them. A 
mechanism which can naturally generate vanishing/small Bfi and inhnite/large 
tan/3 in the MSSM is gauge-mediated SUSY breaking. In its pure form it 
requires Bfi = 0 ai the messenger scale, and a non-vanishing Bfi is only generated 
by renormalization-group running at lower scales. The connection between gauge- 
mediation and the large-tan/3 MSSM has been studied in detail in Ref. 26 


however with several simplifying restrictions which allow for very large but not 
inhnite tan/3. The promising results explained below provide motivation for 
further model building to actually realize inhnite tan /3 in a more fundamental 
theory, but this is beyond the scope of the present paper. 

We now give a brief preview of the MSSM limit tan /3 —)■ cx) to set the stage 
for the remainder of the paper. In the simple case that all relevant SUSY masses 
are equal to Mgusv and tan/3 is moderate, the SUSY contribution to is ap- 
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proximately given by the one-loop diagrams as 


SUSY,1L 


13 X 10 sign(/r) tan (3 


100 GeV\‘ 


Me 


SUSY 


-) 


( 4 ) 


For large tan/3 it has already been shown in Ref. that higher-order terms 
leading in tan /3 can become important and change the linear dependence. These 
higher-order terms can be resummed as 

„SUSY,1L 


,SUSY 


1 + 

where, still in the case of equal SUSY masses, the self-energy leads to 


A cs 


-0.0018 sign(/i) tan/3. 


Hence, in the desired limit of inhnite tan/3 the SUSY contribution becomes 

lim -72 X 10“^° —- , 

tan/3-^oo \MsVSY J 


,SUSY 


( 5 ) 


( 6 ) 


( 7 ) 


demonstrating that the magnitude of the contribution is very large, even for 
Msusy in the multi-TeV range. The proportionality to tan /3 has been replaced by 
a constant behaviour, and the dependence on sign(/i) has disappeared. However, 
the sign is predicted to be wrong, so such a scenario with inhnite tan /3 and equal 
SUSY masses is dehnitely excluded. 

If the SUSY masses are not all equal, the above approximations do not ap¬ 
ply, and in the remainder of the paper we will give a full investigation of the 
hve-dimensional parameter space. In particular we will characterize the regions 
that lead to a positive SUSY contribution to a^, which agrees with the observed 
deviation. 

In Sec.j^we dehne the model, provide analytical results and useful approxima¬ 
tion formulas, and briehy review constraints from other sectors. Comprehensive 
numerical analyses and their physical interpretation for different mass parame¬ 
ter regions are presented in Sec. We also discuss relevant constraints on the 
parameter space. Further discussion and the conclusions are given in Sec. In 
addition, we present relevant benchmark points and scan plots showing the pos¬ 
sible values of the lightest SUSY mass parameter for which the discrepancy ([^ 
can be explained. 


2 The limit tan /3 ^ oo 

2.1 Analytical results and approximations 

The MSSM tree-level muon mass is given by 

= Vi^Vd, 


( 8 ) 
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the product of the muon Yukawa coupling and the down-type Higgs vacuum 
expectation value Vd- Radiative muon mass generation in the MSSM requires the 
tree-level mass to vanish. There are two generic possibilities: either = 0 or 
Ud = 0. In case = 0, the muon mass can be generated from loop diagrams 
involving binos and smuons and the non-standard soft supersymmetry breaking 
The value of in this scenario has already been studied in Refs. 


term vcnuic wi u,^ 

and the observed deviation in a 


20 


21 


^ can indeed be explained for bino and smuon 
masses around 1 TeV. A similar study has been performed with the holomorphic 
trilinear coupling in Ref. 


22 


Here we consider the second possibility = 0, or equivalently the limit of 
Eq. (|^. In the following we explain the relevant formulas governing the muon 
mass, Yukawa coupling, and magnetic moment in this limit. 

The physical muon (pole) mass is given by the on-shell muon self energjj^ 


m, 


= -E 


MSSM _ 




red 

M ‘ 


(9) 


The explicit result can be found in the Appendix. This equation is used to 
determine the value of the muon Yukawa coupling y^. 

Eq. (|^ also dehnes the “reduced” self energy factor (which satishes 


A„ = tan 13 A^f*^ with the usual dehnition for A, 


of Ref. 


24 


and which agrees with 


^). The dehnition highlights the physics behind the 


the quantity 

muon mass generation: chiral symmetry (under which left- and right-handed 
muon transform with different phases) is broken by the Yukawa coupling, and the 
muon mass is generated by loop-induced couplings of the muon to the vacuum 
expectation value of the “wrong” Higgs doublet Vu- Hence the factor y^Vu can be 
pulled out of all contributions to the muon mass. (Another factor of the Higgsino 
mass y could also be pulled out of all contributions, because /i has to appear in all 
couplings of the muon to due to Peccei-Quinn symmetry.) It should be noted 
that the coupling constants, mixing matrices, and mass eigenvalues entering the 
self energy partially depend on the muon Yukawa coupling, so has a small 
residual dependence on y^. 

The combination 




( 10 ) 


is thus an important quantity. In the standard case with moderate tan (3 it is 
strictly equal to tan f3 and also often identihed with tan (3, but the latter 
identihcation is only possible if either higher-order effects can be ignored or on- 
shell renormalization is used. Here we have to distinguish from and 

tan from tan (3. 


^All the following formulas are given for the limit tan/3 —>■ cx). The formulas valid exactly 
also for arbitrary (small, large or infinite) tan /3 can easily be reconstructed by replacing —> 
^ in Eqs. 


13). 
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The physics behind the generation of is the same as for the muon mass. 
To highlight this similarity we write 


^SUSY _ ^red 


( 11 ) 


introducing another dimensionless “reduced” quantity In the standard case, 
we simply have ~ tan/3aj^®'^. Again, has only a small residual depen¬ 

dence on and its explicit result can be found in the Appendix. 

Combining the equations for the muon mass generation and the magnetic 
moment, we obtain the MSSM prediction in the limit tan/3 —)■ oo: 


,SUSY 


,red 


A red ' 


( 12 ) 


We record here further useful relations between the quantities introduced so far: 

1 


,SUSY 


= tan (3 a 


eff ^red 
M ’ 


m, 




tan /3®® = 




A red 


1650//^. (13) 


In contrast to Eqs. 


from the Introduction, the exact result Eq. (12) 


depends on Eve independent SUSY mass parameters: the Higgsino mass y, the 
gaugino (bino and wino) masses Mi 2 , and the left- and right-handed smuon soft 
mass parameters The result has interesting symmetry properties due to 

cancellations between numerator and denominator: on the one hand, is 

invariant under a change of sign of /i and, on the other, it is invariant under 
a simultaneous change of signs of Mi and M 2 . This constitutes a fundamental 
difference compared to the MSSM with finite tan (3, where the sign of can 

directly be specified via the sign of y as e.g. in Eq. (|^. 

Now, useful mass-insertion approximations for the MSSM muon self energy 
and the SUSY contribution to the anomalous magnetic moment of the muon are 
provided. Their validity for the parameter ranges of interest will be discussed at 
the end of the subsequent section. To maximize the magnitude of and minimize 
contraints from CP-violating observables we restrict our considerations to real 
mass parameters p, Mi^ 2 - Taking advantage of the above-mentioned symmetries 
we choose, without loss of generality, /i and Mi to be positive. Only the sign of M 2 
remains arbitrary. Under the assumption Mz Msusy, the SUSY contributions 
are given by the five mass-insertion diagrams in Fig. |^[^[4^. The self energy 
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Figure 1: Mass-insertion diagrams contributing to and For the 

latter, an external photon couples to any of the charged particles in the loop. 


factor can then be decomposed as the sum of 

(14a) 

(14b) 

Af(M/ii)= jLMi/i/(Mi,/x,mL), (14c) 

Af (M/1 h) = - jLMi/i/(Mi,/x,m«), (14d) 

A)f'^(.B/iL/lH) = (14e) 


where gi ^2 are the U(l) and SU(2) gauge couplings, and the loop function is given 
by 


/(a, b, c) 


6^ In |! + 52 c2 In g -f In J 
(a^ — h^) (h^ — c^) (a^ — c^) 


(15) 


This loop function is dimensionful and scales as 1/M^, where M denotes the 
largest of the three mass arguments. The result for the self energy corresponds 
to the one from Ref. 


24 if the limit tan /? —)■ cx) is taken. 


Likewise, for the case Mz <C Msusy, the SUSY contribution to the anomalous 
magnetic moment of the muon can be approximated by the sum of the expressions 
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afiWHD) 

afiWHliL) 

afiBHfiL) 

afiBfiLfiR) 


92 2 

92 2 

' 167r2 

_9]_ 2 

167r2 

9l 2 


Ma/i / M| /i 


m. 


m, 


m 


R 


2 ’ 2 


M 2 /i / M| jj, 


Fh ^ — 

J- h \ 9 5 9 

raf^ 


Mill ^ f Ml Id 


m 2 ^ \^i ’ 

Mi/i ^ fMl 


2 ’ 2 
m^ m|j 


m 2 rn^ 


9l 2 r 

87r2 " V^r Ml 


(16a) 

(16b) 

(16c) 

(16d) 

(16e) 


corresponding to the mass-insertion diagrams in Fig. [TJ with an external photon 
coupling to the charged internal line. Via Eq. (11), this result can be related to 
the one for hnite tan (3 quoted e.g. in Refs. 1^1^, where the loop functions 


Fa{,x,y) 

Fb{x,y) 


Ff (x) - F^{y) 
3>{x-y) 
Fj^jx) - Fj^jy) 
Q{x-y) 


(17a) 

(17b) 


have been introduced. The functions and can be found in the Appendix. 

We now provide seminumerical approximations which allow to directly read 
off numerical orders of magnitude and signs and which facilitate the phenomeno¬ 
logical discussion. The dimensionless loop functions 


/ = ab I {a, b,c), 

(18a) 

KNix,y) = 2xyFb{x,y), 

(18b) 

Kwix,y) = 2xy[2Fa{x,y) - Fb{x,y)], 

(18c) 


are especially useful if the sign and size of the different contributions shall be 
compared. The hrst two of these functions have been introduced in Refs. and 
31 respectively. In the special case of three equal masses, we get 1(1,1) = 1/2, 
Fvi/(1,1) = 5/6 and = 1/6. In general, I and lie between 0 and 

1; Kw lies between 0 and 1.155. Furthermore, in the case of equal arguments, 
I{x,x) and Kn{x,x) increase monotonically, while K\y{x,x) increases up to a 
value of about 1.155 at x ~ 11.15 and then decreases towards a limit of 1 for 

X —>■ oo. 

By inserting all numerically known Standard Model quantities we obtain a 
seminumerical version of Eqs. (14) and the negative MSSM muon self energy. 
















which should be equal to the muon mass [see Eq. (|^], can be written as 

~ I/m 


IM2I 

/i 

mi 

’ rriL^ 

Ml 

jl] 

mx’ 

mi) 

Ml 

j^\ 

rriR 

mRj 

Ml 


rriR 

mn) 


( 19 ) 


where the order of terms is the same as in Eqs. (14), and the first term combines 
the contributions from loops with charged or neutral winos. 

Similarly, by plugging in numbers into Eqs. (16), we obtain the seminumerical 
mass-insertion approximation for in the tan /3 —)■ cx) limit as 


,SUSY 


y^iVu ^red 

(Jj ,, 


m. 


y^l 


sign(M 2 ) 2.483 x 10 


-8 (1 TeV)^ 


0.712 X 10 


1.425 X 10 


1.425 X 10 


/riMal 
-8 ( 1 TeV)^ 

/iMi 

-8 ( 1 TeV)^ 

fiMi 

_8 /xMi (1 TeV)' 


K 


N 


M| 


/i" 


Kw —y, — 
mi mi 

A N I —ir, —^ 

^ ^ 


9 2 


K 


N 


mr m 


R 


Ml' Ml 


( 20 ) 


which corresponds to the result from Ref. in the limit tan/3 —)■ oo and 
combines the contributions that are created by a charged or neutral internal 
wino in the first term. 

The approximate expressions (19) and (20) have several noteworthy features. 

• They are linear in because the residual dependence of and on 

vanishes for M^/Msusy —t 0. 

• All terms involve the factor /iMi or p.M 2 . This explains the symmetry of 

(given by Eq. (12)) under a sign change either of fi or of Mi and M 2 . 

• Under sign change of M 2 alone, the wino contributions change their signs 
relative to the bino contributions. 
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• The signs of the contributions are related: each mass-insertion diagram con¬ 
tributes with equal sign to and except for the chargino diagram, 
where the signs are opposite. 


• Thus if the chargino diagram Q in Fig. [^dominates in both and 
we get <0. In contrast, for neutralino dominance becomes 

positive and, therefore, has the correct sign with regard to an explanation 
of the discrepancy Q. 


The contributions ( |14b[ ) and (16b) from the diagram with an internal neu¬ 
tral wino can never dominate since they are always smaller than the ones 
from the diagram with an internal charged wino. In the charged and 


neutral wino contribution add constructively, in destructively. 


2.2 Experimental constraints from other sectors 


Before we present our numerical analyses, we briefly discuss constraints on pa¬ 
rameter space from other observables and sectors and show that the arising con¬ 
straints can be satished without restricting the hve parameters p. Mi 2 , rel¬ 
evant for a^. Clearly, an observable strongly related to is the lepton flavour 
violating decay p —)■ ey. The correlation of the two MSSM predictions has been 
discussed extensively in the literature 50 -54 . If the prediction for is hxed. 


the SUSY prediction for /i —)■ ey can be estimated quite well up to the unknown 
lepton flavour violating parameters. In particular the one-loop amplitudes for 
the two quantities share the same tan/? enhancement; the correlation between 
the two does not strongly depend on tan/?. E.g. Figs. 13-16 of Ref. 


53 


remain 


valid also for tan/? —)■ 00 . Therefore, like for moderate tan/?, the decay /i —)■ ey 
is compatible with experimental bounds 55 for sufficiently small lepton flavour 
violating parameters. 

Further constraints from i?-physics have been discussed extensively in 
Ref. (26|; taking the limits of those results as tan/? —)■ 00 shows that agree¬ 


ment between theory and experiment can be achieved. The parameters mainly 
constrained by R-physics are the heavy Higgs boson masses as well as the stop 
trilinear coupling At. Since these parameters are not relevant for the discussion 
of a^, we give only a brief account of how they are constrained. The new physics 
contribution to BR(i?+ —)■ r+z/) can be suppressed below 2a with a being the 
experimental error, by having the charged Higgs mass Mh+ around a few TeV 
or higher. We can relax the constraints from Bg —)■ and B —)■ X^y by 

assuming that At and flavour-violating squark mass insertions are small enough. 
In addition, the former process can be further suppressed by raising the masses of 
the UP-odd as well as the CP-even heavier Higgses, and the latter by increasing 
the squark masses. 

Given large bottom and tan Yukawa couplings, large sbottom and stau cor¬ 
rections to the mass of the lighter CP-even Higgs h might undermine the stop 
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contribution that should lift up to the measured value around 125 GeV 56 
as tan/3 —)■ oo. Indeed, Fig. 2 of Ref. shows a rapid drop of nih as tan/3 


grows over around 40. This happens due to a cancellation between the tree-level 
and the contributions to the bottom quark pole mass causing a blow-up of the 
bottom Yukawa coupling. As mentioned in the same reference, this is however 
not necessarily the case for a high tan /3. First of all, it never happens for /i < 0. 
Even if /i > 0, once we enter the regime where | A^l S> 1, the bottom Yukawa cou¬ 
pling comes back to 0 ( 1 ), thereby suppressing its negative contribution to rrih. 
Therefore the tan /3 —)■ cx) scenario allows to choose third generation sfermion 
masses which lead to a correct Higgs mass. 

Owing to decoupling of the multi-TeV extra Higgs masses required by the 
R-physics constraints, the MSSM Higgs sector is approximately SM-like 58 


Nonetheless it is instructive to qualitatively discuss the most relevant SM-like 
Higgs decay modes, h —)■ 66 , h —)■ , and h —)■ which might potentially 

be affected by large bottom and charged lepton Yukawa couplings if the heavy 
Higgs masses are not sufficiently high to apply the decoupling limit. Their decay 
rates can be kept close to the SM values in the following way. The effective 
coupling, for f = b,T, /i, is generated by the same loop diagrams that 
generate my. Therefore, the h-f-f coupling should be SM-like if h consists purely 
of to a sufficient extent. Remember that the coupling, yj, can be 

comparable to or even larger than yt. The remaining question is then how to 
maintain the Higgs mixing angle a small enough for h to avoid excessive H^- 
contamination. It is easy to see that the 2x2 tree-level GP-even Higgs mass 
matrix becomes diagonal as tan/3 —)■ oo (see e.g. Eq. (23) of Ref. [^). The one- 
loop corrections to its off-diagonal elements are suppressed for small At and Ah, 
which can follow from gauge-mediated SUSY breaking at the same time as the 
small By, needed for large tan/3. Regarding h —)■ 77 , we can use the fact that 
the stau-loop contribution (see e.g. Refs. 60-^) to the dimension-hve effective 


operator decouples faster than the tau self energy, as the stau masses increase. 


3 Numerical analysis 

3.1 Dominance transition and sign change 

In the Introduction we have shown that if all SUSY masses are equal, 
becomes negative in the limit of tan /3 —)■ 00 . In this case the wino contributions 


dominate in both and A^®'^, i.e. in the numerator and denominator of Eq. (12), 


and they have opposite signs, see Eqs. (19,20). However, these equations also 


show that the resulting will be positive if any of the three bino contributions 
dominates. 

In the following we will delineate regions in the Eve-dimensional parameter 


space of the masses Mi 2 , y, in which is positive. We start by not- 
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= m|_ = rrip, IMjl = |i = IMgl = p, = rrip 



(a) large /r-limit 


(b) /U/j-dominance 


Figure 2: Individual behaviour of the contributions to and for two dif¬ 
ferent mass spectra. The full result is given by Eq. (12) as the ratio of the 

s um s of the contributions. WHL± denotes the sum oi WHD and WHfiL, for 
positive or negative M 2 , respectively. In the upper plots, the individual contri¬ 
butions to are rescaled by to make the results depend only on ratios 

of SUSY masses. 


ing that the signs only depend on mass ratios and that each of the three bino 
contributions can be expected to dominate if a particular mass hierarchy is valid: 


• The BulUr contributions dominate for /i: “large /i-limit” 


• The BHjjiR contributions dominate for rriL- “/i_R-dominance” 


• The BHfiL contributions dominate for |M 2 |,mij. These have 

the smallest numerical prefactors, and we can expect them to dominate 
only for rather extreme hierarchies]^ 


The phenomenological behaviour and importance of the BfiijiR and BHfiR con¬ 
tributions have also been discussed recently in various contexts in Refs. 31-p5 
and Refs. 34,35,53,63 , respectively. We hrst focus on cases where these 


42,53 


two contributions become dominant. Parameter regions in which the BHJIr con¬ 
tributions dominate will be presented later on. 

Fig. [^illustrates plots of the individual contributions to and A^®*^ for pa¬ 
rameters which interpolate between the equal SUSY mass case and the two mass 


^Dominance of the BHrl contributions can also be achieved with IM 2 I ^ rriL ^ mR. 

We will not investigate this below since the mass hierarchies would be even more extreme. 
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hierarchies for either “large /i-limit”, Fig. [2]Q, or “/i/{-dominance”, Fig. [2]Q. In 
each hierarchy case the desired bino contributions become the largest ones for 
sufficiently large mass ratios. We first discuss Fig. [2]Q, corresponding to the 
“large p-limit”, in detail. It shows the contributions as functions of the ratio 
IM 2 I = fi over Ml = mi = thr] the horizontal axis is given by a; = logiQ^^/niR). 
We pay particular attention to the dominance transition among the contributions 
and the sign change of and 

• Around a; = 0, where all masses are equal, the wino contribution dominates 

in both and For larger values of x, the BjlijiR contributions 

increase proportionally to /i whereas the other contributions are suppressed 
as l//i. Hence the latter tend to zero for large x, except for the WHL con¬ 
tribution to which tends to a constant because it is also proportional 
to M 2 and I M 21 = /i in this plot. 

• At X ~ 0.36, the dominant contribution to changes from the wino 

contribution to the BjlijlR contribution. For larger mass ratios, x ~ 1.00, 
the same dominance change also happens in AJ^®'^ and thus is always 

positive for X > 1. 

• In the intermediate region, 0.36 < x < 1.00, BjlijiR already dominates in 

the numerator but the wino contribution still dominates in the denominator 
of Here the sign of depends on M 2 in such a way that it is 

equal to sign (—M 2 ). 

• For X < 0, the wino contributions remain dominant. 

Similarly, we discuss Fig. [2]Q, which corresponds to “/i/j-dominance”. It 
shows the contributions as functions of the ratio mi over Mi = IM 2 I = p = 
the horizontal axis is given by x = logiQ(m/,/m/{). 

• Like in Fig. [2]Q, the equal mass case emerges at x = 0, and around this 

point the wino contributions dominate in both the numerator and denom¬ 
inator of For larger x, all contributions tend to zero except for the 

BHfiR ones, which are independent of rriR and thus stay constant. 

• At X ~ 0.7, the BHfiR contributions become dominant; accidentally, the 

dominance change occurs approximately at the same x in the numerator 
and denominator. For larger values of x, this leads to a positive , 

independently of the sign of M 2 . 

To investigate the sign of comprehensively we plot it in a more general 

parameter space, as a function of the two mass ratios used in Figs. |2]Q and 
@, he. in the plane rriL/mR versus n/mR. Fig. [^displays the results for four 
different choices of the remaining mass ratios Mi/ruR and |M 2 |/p. We begin by 
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Figure 3: Signs of in the plane of versus fi/ruR, for the two possible 

signs of M 2 . In the white regions is positive. The red regions indicate 

negative for M 2 > 0, the blue hatched ones negative for M 2 < 0. 


In the overlap regions 


,SUSY 


is always negative. The first plot also displays 


which contributions dominate in the bulk of each region in the numerator and 
denominator of Eq. (12). (It should be noted that near the borders between the 
regions there are cancellations, and different contributions can be the largest.) 
In each plot the remaining two mass ratios Mi/rriR and IM 2 I//U are fixed as 
indicated; the structure of the dominance regions is the same in all plots. For the 
small white region in the bottom left area of plot (d), see the text. 
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discussing Fig. [sjQ. It generalizes Fig. |^by assuming Mi = uir and IM 2 I = /i, 
so the different regions in Fig. [sjQ can be fully understood from the previous 
discussion. 

• The large, central red/blue overlap region contains the origin at which all 

relevant SUSY masses are equal. In this region the wino contributions 
dominate and is always negative. 

• The two white regions, where is always positive, correspond to the 

“large /i-limit” and the “pR-dominance” region. 

• The blue hatched and the red regions (excluding the red/blue overlap) are 

transition regions where different contributions dominate in numerator and 
denominator. They generalize the region 0.36 < x < 1.00 of Fig. [2]Q. Here 
the sign of depends on M 2 and is equal to the sign of either (+M 2 ) 

(blue hatched) or {—M 2 ) (red). 

In Fig. |3]Q, the gaugino masses IM 2 I and Mi are equated with p and 
respectively. This is a reasonable choice since their precise values are qualitatively 
unimportant for the appearance of the different regions. The dependence on 
these two gaugino masses is investigated in Figs. ag-Q. whose presentation is 
similar to Fig. [sjQ but with either Mi = 10IM 2 I = 10/i, or Mi = 0.1 
and IM 2 I = 100/i. Together, these plots cover large, representative parts of the 
full parameter space of the four relevant mass ratios. 

Obviously, all plots can be understood in the same way as Fig. |3]Q. As Mi 
becomes heavier than in Fig. |3]Q, the bino contributions are suppressed, the 
wino-dominance area expands, and the positive regions shrink. In contrast, 

the positive regions expand in Figs. |3]0 and (j^, since M 2 gets heavier 

than /i and the wino contributions are thereby suppressed. 

Finally we investigate the positive region with BHfii dominance, which 

cannot be seen for the most part in the previous graphs. This dominance can be 
realized under the hierarchy condition Mi, p, rriL IM 2 I, which is satisfied in 
the small white region in the bottom left area of Fig. [sjQ. To better understand 
the BHJIl dominance, we change the plot axes to the ones of Fig. ijQ. Here the 
sign of as well as information about which contributions dominate in 

and are displayed in the plane of mu/^ versus Mi//i. Fig. 4 b| corresponds 
to the y-axis of Fig. |4jQ and shows the individual behaviour of eacn contribution. 
As we can see, the BH/ll contributions dominate only for extreme mass ratios 

mfl/zi > 100. 


3.2 Magnitude of and constraints on Yukawa cou¬ 


plings 


Having understood the parameter regions which lead to positive we now 




turn our interest to its magnitude and also take into account other constraints 
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Figure 4: Q Signs of in the plane of m^/fi versus Mi//x, for the two 

possible signs of M 2 . The colour coding is the same as in Fig. (|^ Individual 
behaviour of the contributions to and as a function of the mass ratio 
mn/jj,. The presentation is analogous to Fig. |2| 


on the parameter space. As discussed in the Introduction, we can write 

ar > ^^ = (^(1), (21) 

SUSY,min 

where MsusY.min is the smallest among the hve relevant mass parameters Mi, 
I M 2 1, /i, and mL,R [and thus similar but not equal to the lightest SUSY particle 
(LSP) mass since the latter arises from diagonalization of mass matrices]. The 
coefficient C is expected to be of order unity because of radiative muon mass 
generation. It depends only on the four independent mass ratios (up to terms 
suppressed by powers of M|/M|ugY)- Hence for each given set of mass ratios we 
can pose two questions: 


1. What is the value of C? 

2. What is the minimum SUSY mass for which = 28.7 x 10“^°? 


The answers are given in Fig. [^for the special choices Mi = thr and IM 2 I = fi 
which have already been used in Fig. [sjQ. The hgures also show the behaviour 
of the muon Yukawa coupling, which is proportional to see Eq. (13). The 

four most interesting regions are the “large p-limit” and “/iij-dominance” regions 
with either positive or negative M 2 . These are located in the right end and middle 
top areas of the two plots and correspond to the white regions of Fig. |3]Q with 
a positive In detail Fig. shows the following. 
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I'^SUSY.min [TeV] 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



0.62 0.88 1.08 1.25 1.39 1.53 1.65 1.76 1.87 1.97 



Figure 5: Magnitude of and muon Yukawa coupling in the plane of 

versus y/mn, for positive M 2 (left) and negative M 2 (right). The colour 
coding corresponds to the different values of the coefficient C defined in Eq. (21) 
and the equivalent minimum SUSY mass which leads to agreement with the 
current deviation via = 28.7 x 10“^°. The black contours correspond to the 

indicated values of the muon Yukawa coupling resulting from Eq. (g. The grey 
regions are excluded by negative In the white regions the denominator 

undergoes a sign change, and the perturbation theory becomes untrustworthy due 
to large Benchmark parameter points introduced in Table are marked with 
blue squares. 


• The first important observation is that in the centre of the “large /i-limit” 
and “/iij-dominance” regions, C lies between 0.2 and 0.4; the correspond¬ 
ing minimum SUSY mass is around 1 TeV. C does not quite reach unity 
because of the behaviour of the loop functions in and 

• In both regions, the magnitude of U, the minimum SUSY mass, and the 
muon Yukawa coupling depend on the sign of M 2 . This dependence arises 
from interference of the dominant contributions with subdominant wino 
contributions, which can be either constructive in A))®'^ and destructive in 
^red QJ, versa. 

• There are small white strips in parameter space with larger C and minimum 
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Figure 6: Excluded parameter regions in the plane of mi/m^ versus 
for positive M 2 (left) and negative M 2 (right). The areas corresponding to the 
different constraints are plotted in the same order as they appear in the legend. 
Detailed explanations for each constraint can be found in the text. Benchmark 
parameter points introduced in Table are marked with blue squares. 


SUSY masses. In these strips undergoes a sign change and the Yukawa 
coupling becomes inhnite. Here perturbation theory is not trustworthy since 
two-loop effects are non-negligible due to the accidental cancellation of the 
one-loop contributions. We exclude these regions from further discussion. 


Of particular interest are the regions with smaller Yukawa couplings, i.e. 
the “large /i-limit” region for negative M 2 and the “/i/j-dominance” region 
for positive M 2 . Here C can still increase up to 0.3, and the minimum 
SUSY mass reaches up to 1.1 TeV. The opposite sign choices lead to the 
strips where diverges. In the “large /i-limit” the muon Yukawa coupling 
is always proportional to l//i and therefore its magnitude is comparatively 
small and decreasing as y increases. In the regions of “/lij-dominance” \y^\ is 
generally much larger. Eq. (19) allows to deduce that for “/i/j-dominance”, 
ly^l cannot become smaller than 0.78 for negative M 2 and 0.12 for positive 
M 2 , respectively. 
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Now, we consider constraints on the relevant parameter space for large 
All previous plots are based on the mass-insertion approximation and depend 
only on mass ratios. For each given set of mass ratios with positive there is 

a value of the overall mass scale Mgusv for which agrees with experiment. 

In the following discussions we always £x Mgusv to that particular scale. 

Fig. § shows the same parameter space as Fig. but displays regions in 
which charged SUSY masses become too small (medium grey), the muon Yukawa 
coupling is non-perturbative (dark grey) or the electroweak vacuum becomes 
unstable (hatched). In detail, the constraints are the following: 

• The light grey regions, which have already been shown in Fig. yield 

negative values of and are thus not of interest with regard to an 

explanation of the discrepancy ([^. 

• There are strong but model-dependent collider bounds on the masses of 
charged SUSY particles. Here we are interested in very large electroweak 
SUSY masses, which are far above the LHC limits. Therefore we show as 
the medium grey regions in Fig. [^an exemplary contour corresponding to 
a chargino or smuon which is lighter than 100 GeV. 


In the dark grey regions the muon Yukawa coupling violates perturbativity, 
\y^\ > \/47r. Generally, we regard the model as a pure low-energy model 
and do not require the stronger constraint that there be no Landau poles at 
higher scales; such an analysis has been done in Refs. 1^,26 


It should be 

further noted that the observable couplings of the lightest Higgs boson are 
SM-like up to corrections suppressed by powers of the SUSY scale. This is 
true independently of the value of the fundamental Yukawa coupling due 
to decoupling. 


As discussed in Refs. |32[|60f|62|64| , the combination j/^/i is limited by the re¬ 
quirement that the electroweak vacuum is metastable with a lifetime longer 
than the age of the universe. This constrains in particular the parameter 
regions with large /i as already stressed in Ref. (^. The regions which are 
excluded by the htting formula from this reference applied to the muon sec¬ 
tor are shown in Fig. Since it is possible that yi > yt in our scenario, the 
0{y^) counterpart of the standard 0{y^) correction to the effective poten¬ 
tial might be signihcant. To scrutinize the effect of the former, we add to the 
effective potential the one-loop contributions (see e.g. Eq. (1.2) of Ref. [65] ) 
in the gaugeless limit, with the renormalization scale Q = They 

include the (s)lepton- and the Higgs(ino)- as well as the (s)top-loop con¬ 
tributions. We then evaluate the Euclidean action numerically using the 
method presented in Ref. 66 . In principle, the metastability bound on 


yifi depends on the pseudo-scalar Higgs mass Ma affecting the one-loop 
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Table 1: Benchmark parameter points. All masses are given in TeV. 


correction to the slepton squared mass. Motivated by the i?-physics con¬ 
straints, we set Ma = 50 TeV, a high enough value for extra Higgs states 
to decouple. By applying the specihed procedure to the muon sector, we 
hnd that the regions shown in Fig. [^are excluded because of an unstable 
electroweak vacuum. 


According to the results shown in Fig.|^ the requirement of vacuum metastability 
poses strong constraints on parameter space, in particular on the “large /i-limit” 
region. This region is completely excluded for positive M 2 , while for negative 
M 2 the ratio ii/rrtR must not be too large. The dominance” region is not 
constrained by vacuum metastability. The two different approaches to the vac¬ 
uum metastability yield similar results but the influence of the additional terms 
is visible in Fig. [6]Q. 

The Yukawa coupling of the r lepton, which is another quantity of interest, 
can be determined similarly to y^. There is a strong correlation if the stau 
masses are assumed to be equal to the smuon masses; without this assumption 
the two Yukawa couplings are basically independent, and vacuum metastability 
from the tan sector does not further constrain the hve-dimensional parameter 
space. However, a detailed analysis should take into account higher orders in the 
large tan Yukawa coupling and is beyond the scope of this paper. 

Thus we hnd that out of the four interesting regions in Fig. three survive. 
Now we briehy comment on the comparison between the approximate results for 
as well as obtained from Eqs. (14,16) and their exact results according to 
the Appendix. In the three viable regions described above, the deviation amounts 
to at most a few percent. Outside these regions the approximation becomes worse 
for very small SUSY masses and in the small strips with divergent muon Yukawa 
coupling. 


4 Conclusions 

We have studied the muon magnetic moment in the MSSM for tan/? —)■ 00 . This 
is a viable limit in which all down-type masses are generated by loop-induced 
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couplings to the “wrong” Higgs Hu- The scenario can also be viewed indepen¬ 
dently of supersymmetry, as a simplihed model which realizes the general idea of 
radiative muon mass generation. 

The SUSY contribution to the muon anomalous magnetic moment is given by 
Eq. (12) as the ratio of two one-loop quantities. For this reason the scenario leads 
to large with TeV-scale new physics masses, and it behaves qualitatively very 
differently from the standard MSSM with tan /3 < 50. E.g. the signs of gaugino 
and Higgsino masses drop out between numerator and denominator and cannot 
be used to adjust the overall sign of For equal masses, is strictly 

negative, see Eq. Q. In general, the sign of depends on the ratios between 

the hve relevant mass parameters; Figs. [2]-|^give a comprehensive explanation of 
which mass ratios lead to positive 

The three most promising parameter regions with positive are the 

“large /i-limit” for negative M 2 and the “/f/j-dominance” region for positive or 
negative M 2 . The “large /x-limit” region with positive M 2 is excluded because 
of vacuum instability. In all three viable cases, large contributions to are 

possible without conflicts with other experimental constraints, see Figs. The 
discrepancy Q can be fully explained with a lightest SUSY mass above 1 TeV. 
If this scenario is correct, the LHC would not be able to hnd the SUSY particles 
relevant for a^. However, Sec. 2.2 shows that the scenario could potentially in¬ 
fluence and thus be tested in a multitude of observables in H-physics and Higgs 
physics. 

To conclude, we summarize the behaviour in three simple ways. Table lists 
several benchmark parameter points with SUSY mass parameters at 1 TeV or 
higher and gives the resulting and muon Yukawa coupling. Two points cor¬ 

respond to the “large /x-limit” with negative M 2 and three to the “/iij-dominance” 
region with positive or negative M 2 . Some points involve only one heavy and four 
lighter SUSY masses; for these large is possible for less extreme mass split¬ 

ting than for the ones with two heavy and three lighter masses. The benchmark 
points further illustrate that the Yukawa coupling is much larger in the case of 
“yhij-dominance”. 

Fig.|7| shows the results of two scans of the hve-dimensional parameter space: 
plotted are the largest possible values of the minimum SUSY mass MsusY,min 
which lead to agreement with Eq. Q . The left panel shows the results as a func¬ 
tion of log^o(^i/^i) positive M 2 ), it corresponds to the y-axis of Fig. |5]Q; 
the right panel as a function of logio(At/Mi) (for negative M 2 ), it corresponds to 
the x-axis of Fig. |5]Q. The results correspond to the regions presented in Fig. 
and are just slightly higher than the values obtained there with hxed mass ratios. 

Finally, for the regions of interest, a simple approximation in the style of 
Eq. ([^ can be derived: for the “large /i-limit” region, we consider the limit 
|/i| S> |Mi| = rriL = mR = Msusy, and for the “/iij-dominance” region we 
consider the limit mi S> |/i| = |Mi| = uir = Msusy- In both of these limits, we 
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Figure 7: Results of scans over all five relevant SUSY mass parameters. The 
plots show the maximum values of MgusY,min for which the deviation Q can 
be explained, also indicating the characteristic of the parameter point for which 
the maximum is achieved. In some regions the achievable maximum depends 
on the cutoff on the magnitude of the muon Yukawa coupling, as indicated; in 
the other regions the maximum is achieved for smaller Yukawa couplings. The 
scans take into account the constraint of vacuum metastability (one-loop), as 


discussed in Sec. |3.2[ and in the right plot all investigated parameter points with 
XogiQ^n/Ml) > 1.8 are excluded by this. 
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A Appendix 


In this Appendix we provide the full one-loop results for the muon self energy and 
the magnetic moment, defined in Eqs. ([9, 11). We use the notation of Ref. 35 


which is similar to the ones from Refs. 28,^. Note that in the standard case 
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of the MSSM with moderate tan/3 < 50, the combination m^tan/3 frequently 
appears. An important difference is that here we have to take into account the 
distinction between the tree-level and the physical muon mass, 7 ^ m^, and 
we have to use the identihcation ( 10 ), 


tan (3 = tan = y^i 


v/2M, 


w 


92 


(23) 


which is exactly valid. Here (72 denotes the SU(2) gauge coupling. 

The relevant interaction terms for chargino-sneutrino and neutralino-smuon 
interactions with muons can be written as 


= Xi + c^u^Pr) + {n%^PL + nf-^^PR) /x/ll + h. c., (24) 

where the coupling coefficients are dehned as 


= -92Va^ 

= y^^Ui2, 

A. = (ffiXi + S2Nk) K - V.NhUL 


(25a) 

(25b) 

(25c) 

(25d) 


The unitary matrices U, V, N diagonalize the chargino and neutralino mass ma¬ 
trices; is used to diagonalize the smuon mass matrix. The resulting smuon 
mass eigenvalues are given by 


= 2\^L + 'mR + 


Ml 


± 


m, 


m 


R 


+ Ml - 


2 e 2 


1 2 


+ 8yl\ 


\2 ^W 

9 I 


(26) 


A term containing the trilinear coupling vanishes because of Vd = 0, and in 
the last term we have used Eq. (23). The occurring SUSY parameters are the 
higgsino mass parameter p, the gaugino masses Mi as well as M 2 and the soft 
mass parameters for the second-generation slepton doublet and singlet and 
rriR, respectively. 

With this notation, the MSSM one-loop contributions to the on-shell muon 
self-energy from chargino-sneutrino or neutralino-smuon loops can be expressed 
as Sf SM ^ sf + Sf = with 


_ _ 

167r2 


sf = 






IGvr^ 


^m^oRe 


(27a) 

(27b) 


j,k 
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where terms suppressed by m^/MsvsY or by powers of tan/? are neglected. In 
dimensional regularization with parameter e = {4 — D)/2 and regularization scale 
fiR, the one-loop integral Bq is given by 


Boi^a, 6) — A -|- 1 -|- 
with the common quantity 


1 


6 ^ — 


a? In 




+ 


(28) 


A = —h In dvr — 7 e, 
e 


(29) 


containing the Euler-Mascheroni constant 'Je- Owing to the unitarity of the 
mixing matrices, in Eqs. (27) the UV divergences as well as the dependence on 
vanish. 

Because of Vd = 0, the tree-level mass of the muon (as well as the other 
charged leptons and down-type quarks) equals zero. Therefore, the entire muon 
mass has to be generated at the loop level, which means 




MSSM 


(30) 


For a given parameter point, this equation can be used to determine the muon 
Yukawa coupling Note that the self energy is approximately linear in y^, as 
suggested by the dehnition of A)f'^ in Eq. (|^. But the Yukawa coupling also 
enters nonlinearly, through quadratic terms in n^- {nfr, )*, and via the smuon 
masses (26). Therefore, it has to be determined numerically in an exact calcu¬ 
lation. As shown in Sec. however, the nonlinear terms vanish in the mass- 
insertion approximation. 

The SUSY one-loop diagrams contributing to the anomalous magnetic mo¬ 
ment of the muon are the same as for the muon self-energy, except for the fact 
that now a photon couples to the charged SUSY particle in the loop, i.e. the 
chargino or the smuon. Again neglecting terms suppressed by tan/?, the result 

= ( 2 /;,u„/m^)ayf'^, with 


reads = a 


f+«? 


af = 


at = 


247r^ 


m, 
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9 

mi 


Re 


c^- 


(cSj- 


Cl 


m^o 

A 




dSTT^ mi 

j,k 

where the mass ratios are given by Xi = m‘i_/mi and xjk = mio/m 
dimensionless loop functions have been introduced in Ref. 28 as 

3 


2 

i^k 


Fy 


(a:) = 


?N, 


X) = 


2(1 -x)3 
3 

(1 -a:)= 


[—3 -I- 4a: — — 2 In x] , 

[l — a:^ -I- 2a:lnx] , 


(31a) 

(31b) 

The 

(32a) 

(32b) 
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normalized such that -F^(l) = -^2^(1) = 1- The formulas look identical to the 
respective terms in the standard case with moderate tan /d, but like for the self 
energy, there are nonlinear terms in which usually are neglected because 
of the m^/MsusY suppression. In our case, these terms are only suppressed 
as Mw,z /hence they are much more signihcant in an exact calculation. 
Nevertheless, like for the self energy, in the mass-insertion approximation, 
is linear in the Yukawa coupling. 
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